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ABSTRACT 

A new technique has been developed to remove the ionosphere's distorting effects from low 
frequency VLBI data. By fitting dispersive and non-dispersive components to the phases of multi- 
frequency data, the ionosphere can be effectively removed from the data without the use of a 
priori calibration information. This technique, along with the new gating capability of the VLBA 
correlator, was used to perform accurate astrometry on pulsar B0950+08, resulting in a much 
improved measurement of this pulsar's proper motion = —1.6 ± 0.4 mas/yr, = 29.5 ± 0.5 
mas/yr) and parallax (tt = 3.6 ±0.3 mas). This puts the pulsar at a distance of 280 ±25 parsecs, 
about twice as far as previous estimates, but in good agreement with models of the electron 
density in the local bubble. 

Subject headings: astrometry — techniques: interferometric — pulsars: individual (PSRB0950+08) 



1. Introduction 

Accurate radio pulsar distances are important 
for population and birthrate modeling and, com- 
bined with measurements of dispersion and Fara- 
day rotation, for studies of the ionized interstellar 
medium and Galactic magnetic fields. Parallax 
measurements can firmly fix the bottom rungs of 
the pulsar distance ladder, but parallaxes of even 
the closest pulsars are only a few milliarcseconds 
in amplitude, and such measurements are avail- 
able for only nine pulsars (see Toscano et al. 1999 
for a recent list). 5 

Four of these parallax measurements were made 
using timing studies of millisecond pulsars with 
particularly stable spin periods. The others were 
made with various VLBI techniques. The use of 
calibration sources within the primary antenna 
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An optical parallax is also available for the pulsar 
B0633+17 (Geminga), which has no known radio counter- 
part. 



beam shows promise, as demonstrated by Foma- 
lont et al. (1999). Observing at L band with 25 
meter telescopes provides a primary beam size of 
about 30 minutes of arc. Only a small fraction of 
pulsars have calibrators within this small distance. 
The number of accessible target pulsars is greatly 
increased by the use of calibrators at greater off- 
set angles, up to several degrees, but differential 
ionospheric phase delays make high precision com- 
parison measurements very difficult. 

Most pulsars are weak radio sources at high fre- 
quencies, where array resolution is high and iono- 
spheric effects are minimized. Observations are of- 
ten made at L band (1.4 to 1.7 GHz) as a compro- 
mise, but even at L band the ionosphere remains 
a significant obstacle. We have developed a new 
calibration technique to remove the effects of the 
ionosphere from the data and used it to measure 
a revised parallax of PSR B0950+08. 

2. Observations 

At L band, the VLBA provides simultaneous 
access to a large spanned bandwidth and is sensi- 
tive enough to detect a pulsar within the changing 
ionosphere's timescale of about one minute. Re- 
cently, gating the VLBA correlator became pos- 
sible. This allows the correlator to be turned off 
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during the off pulse phase of the pulsar, increasing 
the signal to noise of the pulsar by a factor of 3 
to 6, depending on the pulse duty cycle. For PSR 
B0950+08 the improvement was a factor of 3.6. 

We observed the pulsar simultaneously in eight 
8 MHz bands ranging from 1410 MHz to 1730 
MHz. Each band was split into thirty two 250 
kHz channels. The VLBA antennas nod back and 
forth between the pulsar and the calibrator with 
cycle times of about five minutes. Seven epochs 
from 1994 to 1999 were observed near parallax ex- 
trema. The first three epochs were ungated and 
used 1004+141 as a calibrator source. Its dis- 
tance of seven degrees from the pulsar made phase 
referencing difficult. These data were useful for 
developing techniques and understanding limita- 
tions but were not of high enough quality for ac- 
curate astrometry. The more recent epochs used 
J0946+1017, less than three degrees distant from 
PSR B0950+08, as a calibrator. 

Most astrometry at low frequency is done at 
night since the ionosphere phase fluctuations are 
two or three times stronger during the day. Since 
parallax measurements are most effective when the 
target and the sun are 90 degrees apart in the 
sky, these observations need to be done at sun- 
rise and sunset — the times when the ionosphere 
changes most rapidly. 

3. Data Reduction 
3.1. Phase-referencing 

This is a VLBI phase-referenced (Taylor et al. 
1999) experiment with an additional calibration 
step. A strong source, J0946+1017, was used as 
the reference source. Using this object's known 
position, the fringe rates and delays were mea- 
sured and interpolated over the pulsar scans. The 
calibrator is self-calibrated so that its structure 
can be modeled. The self-calibration phases are 
interpolated over the pulsar. This simple phase 
referencing determines the pulsar's position to less 
than 30 mas. The phase center is shifted to 
the pulsar's position and the 32 spectral channels 
within each 8 MHz band are averaged together. 
The large number of channels prevents bandwidth 
smearing from decorrelating the long baseline vis- 
ibilities. 



3.2. The Ionosphere 

The ionosphere, like the free electrons in the 
ISM, introduces a frequency dependent delay to 
the radio waves. Unlike for the ISM, this delay can 
be drastically different above each antenna since 
the ionosphere has structure on all scales, ranging 
from meter-sized variations to the bulge caused by 
the sun's radiation which stretches a third of the 
way around Earth. The ionosphere is also differ- 
ent along different lines of sight from the same an- 
tenna. Since the lines of sight from each telescope 
to the pulsar and calibrator puncture the iono- 
sphere at different points, phase-referencing only 
calibrates away the bulk ionosphere over each an- 
tenna, leaving spatial gradients in the ionosphere 
uncorrected. These gradients can change signifi- 
cantly on time scales shorter than one minute and 
are large enough to move the pulsar by as much as 
10 mas per degree of source separation. A phase 
referenced image with no ionospheric correction is 
shown in figure 1. 

At a given time on a given baseline, the phase- 
referenced visibility phases have this structure: 
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non-dispersive 



dispersive 



U and V are the projected baseline lengths, I and 
m are the pulsar's vector distance from the phase 
center (the quantity that this experiment is try- 
ing to measure), and 5 A v is the phase due to 
the unmodeled wet troposphere and is very small, 
introducing less than 45 degrees of phase. These 
three terms are non-dispersive and hence have the 
same frequency dependence. 5 A is measured in 
rad/GHz, and SB is in rad-GHz. SB /vis the phase 
due to uncorrected ionosphere. Since the remain- 
ing ionosphere is the only dispersive term, it can 
be separated from the non-dispersive sources of 
phase. 

This separation, however, is complicated by 
phase ambiguity. For a given set of phase measure- 
ments, a whole family of solutions can be formed 
since n-2ir can be added to the original phases. Be- 
cause our bandwidth is limited, any integral value 
of n will generate a different solution with an al- 
most equally good fit. It is fortunate that the value 
of 5 A is always small enough that a good guess of 
the pulsar's position constrains the ambiguity to 
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a single value of n. 

3.3. Determining the Pulsar's Position 

A bootstrapping approach is used to determine 
the pulsar's position. Using just the 5 Southwest- 
ern VLBA antennas ( Los Alamos NM, Pie Town 
NM, Fort Davis TX, Kitt Peak AZ, and Owens 
Valley CA ) with a maximum baseline of 1500 km, 
a separate image is made for each 8 MHz band in 
the phase-referenced data. Figures 2 and 3 show 
two such images. Plotting the measured positions 
as a function of frequency reveals a trend - the 
points line up in frequency order. If we assume 
a simplified model of the ionosphere as a single 
wedge lying above all the antennas then we would 
expect the pulsar to appear at P{v) = P 00 +w/v' 2 . 
Fitting for P^ , the extrapolated infinite frequency 
position, yields an initial guess of the pulsar's po- 
sition ( the star on figures 2 and 3 ) with a pre- 
cision of 3-10 mas. The wedge model is certainly 
far from correct over the whole Earth but works 
fairly well for the Southwest array. Because the 
ionosphere changes so quickly, a 15 to 30 minute 
snapshot usually works better than imaging all of 
the data. 




Fig. 1. — Image of pulsar before removing the 
ionosphere. All of the flux is from the pulsar, scat- 
tered by the ionosphere. The beam is 4 by 11 mas. 
The lowest contour level is at 5 mJy and a factor of 
2 separates contours. All data from all telescopes 
were used to produce this image. 
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Fig. 2. — 30 minute snapshot made at 1.40 GHz 
using only the 5 Southwestern antennas. The 
beam size and shape match that of the second- 
highest contour. The three pluses designate mea- 
sured location of the pulsar at three different fre- 
quencies; from the top, these frequencies are: 1.40 
GHz, 1.63 GHz, and 1.73 GHz. The star indicates 
the inferred true location of the pulsar by extrapo- 
lating the measured positions to infinite frequency. 
The center of the square is where the pulsar was 
determined to be using the techniques shown in 
this paper. 



To get antenna-based phases as a function of 
frequency, the pulsar is self-calibrated with a so- 
lution interval of between twelve and thirty sec- 
onds and the antenna-based phase corrections are 
recorded. The self-calibration phases cannot be 
applied to the pulsar, as that would move the pul- 
sar to the phase center. For each of the South- 
western antennas, the nondispersive and disper- 
sive components of the self-calibration solutions 
are fit, holding the geometric £ and m terms fixed. 
This is reasonable since we know the pulsar's loca- 
tion to much less than 30 mas, the fringe spacing 
of the longest Southwestern baseline. The other 
antennas arc initially ignored as the phase contri- 
bution due to geometric delay is significant on the 
longer baselines. A new calibration table is made 
containing only the dispersive 5B/v corrections 
that were fit to the data. Imaging the data with 
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Fig. 3.— Same as figure 2 but at 1.73 GHz. 

these corrections produces a new position accurate 
to less than a milliarcsecond. With the improved 
pulsar position, the fitting for SB and imaging are 
repeated with all antennas. The final calibration 
table contains phase corrections for each antenna 
tabulated at time intervals of a minute or less. 
Figure 4 shows the SB values for the Kitt Peak- 
to-North Liberty baseline. This final image has 
much less scattered power and all frequencies ap- 
pear at the same location (figure 5). 
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In summary, although the self-calibration pro- 
cess is used to derive antenna-based phases over 
short periods of time, only the dispersive part of 
these phase terms is removed, removing the iono- 
spheric distortions. The non-dispersive part of the 
phases which contain information about the source 
position (as well as small atmospheric terms) re- 
main in the data. 

4. Application of Ionospheric Removal to 
PSR B0950+08 

PSR B0950+08 has a published parallax of 
7.9 ± 0.8 mas (Gwinn 1986). This measurement 
was made with a three element VLBI array con- 
sisting of Arecibo, the 140 ft. antenna at Green 
Bank and the 40 m antenna at Owens Valley. A 
3-element interferometer produces a beam pattern 
with severe lobe ambiguity. To break the degener- 
acy of the solution, the proper motion was used to 
select a self-consistent set of lobes. Unfortunately, 
the published values of /i Q = 15 ± 8 mas/yr and 
Us = 31 ± 4 mas/yr (Lyne et al. 1982) were not 
sufficiently accurate, causing an incorrect choice 
of lobes. 

I i I i i I 
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Fig. 5. — Image of the pulsar after removing the 
ionosphere. The same region as figure 1 is shown 
with the same contour levels. The highest contour 
is at 64 mJy 



Fig. 4. — Time variation of SB on the Kitt Peak 
to North Liberty baseline. 



The four gated epochs of data were calibrated 
as described and fit with a model with five com- 
ponents: position of pulsar (two components), 



proper motion of pulsar (two) and parallax (one). 
The fit is shown in figure 6, and the proper motion 
and parallax are: 

fi a = — 1.6 ± 0.4 mas/yr 
yU5 = 29.5 ± 0.5 mas/yr 
7T = 3.6 ±0.3 mas 

The first three epochs of data suggest a parallax 
of 4 ± 1 mas suggesting some stability of correc- 
tion even with a calibrator separation of 7 degrees. 
The nutation model in the correlator may leave an 
annual signature with an amplitude of about 40 
/zas per degree of separation (C. Ma et al. 1998), 
about 0.1 mas in this experiment. Systematic er- 
rors such as these are not included in the error 
bars but are expected to be less than 0.3 mas. 
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changing the local electron density from n e = 
0.023cm -3 to n e = 0.01cm -3 , a value much closer 
to the estimated electron density inside our local 
bubble, n e rs 0.005cm -3 . 
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Fig. 6. — Location of pulsar at four epochs and 
model fit to the points. 



5. Implications 

There has been a long-standing problem of rec- 
onciling the electron density derived from the dis- 
tance to PSR B0950+08 with that based on X- 

ray data (Toscano et al. 1999). This measure- This 2 - C ol U mn preprint was prepared with the AAS IATeX 

ment puts the pulsar more than twice as far away, macros v5.0. 
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